import numpy as np
import matplotlib.pyplot as plt

A1 = np.loadtxt('100K_25meV.txt')
A2 = np.loadtxt('200K_25meV.txt')
A3 = np.loadtxt('300K_25meV.txt')
A4 = np.loadtxt('400K_25meV.txt')
A5 = np.loadtxt('500K_25meV.txt')
A6 = np.loadtxt('600K_25meV.txt')
A7 = np.loadtxt('680K_25meV.txt')

plt.errorbar(A1[5:,0],A1[5:,1],yerr=A1[5:,2],marker='o',markersize=6,label='100K',color='navy')
plt.errorbar(A2[5:,0],A2[5:,1]+0.05,yerr=A2[5:,2],marker='s',markersize=6,label='200K',color='blue')
plt.errorbar(A3[5:,0],A3[5:,1]+0.10,yerr=A3[5:,2],marker='p',markersize=6,label='300K',color='cyan')
plt.errorbar(A4[5:,0],A4[5:,1]+0.15,yerr=A4[5:,2],marker='v',markersize=6,label='400K',color='lime')
plt.errorbar(A5[5:,0],A5[5:,1]+0.20,yerr=A5[5:,2],marker='d',markersize=6,label='500K',color='gold')
plt.errorbar(A6[5:,0],A6[5:,1]+0.25,yerr=A6[5:,2],marker='h',markersize=6,label='600K',color='red')
plt.errorbar(A7[5:,0],A7[5:,1]+0.30,yerr=A7[5:,2],marker='X',markersize=6,label='680K',color='darkred')

plt.text(0.2,0.01,'100K',color='navy',fontsize=20,ha='left',va='bottom')
plt.text(0.2,0.06,'200K',color='blue',fontsize=20,ha='left',va='bottom')
plt.text(0.2,0.11,'300K',color='cyan',fontsize=20,ha='left',va='bottom')
plt.text(0.2,0.16,'400K',color='lime',fontsize=20,ha='left',va='bottom')
plt.text(0.2,0.21,'500K',color='gold',fontsize=20,ha='left',va='bottom')
plt.text(0.2,0.265,'600K',color='red',fontsize=20,ha='left',va='bottom')
plt.text(0.2,0.32,'680K',color='darkred',fontsize=20,ha='left',va='bottom')

plt.plot(A1[5:,0],A1[5:,1]+0.05,ls='dashed',label='_nolegend_',color='navy',zorder=0,lw=1)
plt.plot(A1[5:,0],A1[5:,1]+0.10,ls='dashed',label='_nolegend_',color='navy',zorder=0,lw=1)
plt.plot(A1[5:,0],A1[5:,1]+0.15,ls='dashed',label='_nolegend_',color='navy',zorder=0,lw=1)
plt.plot(A1[5:,0],A1[5:,1]+0.20,ls='dashed',label='_nolegend_',color='navy',zorder=0,lw=1)
plt.plot(A1[5:,0],A1[5:,1]+0.25,ls='dashed',label='_nolegend_',color='navy',zorder=0,lw=1)
plt.plot(A1[5:,0],A1[5:,1]+0.30,ls='dashed',label='_nolegend_',color='navy',zorder=0,lw=1)


plt.minorticks_on()
plt.xlim(0,20)
plt.ylim(0,0.44)
#plt.legend(frameon=False,fontsize=14)
plt.xlabel('Energy (meV)',fontsize=20)
plt.ylabel('DOS (1/meV)',fontsize=20)
plt.tick_params(which='both',direction='in',labelsize=20)

plt.savefig('25meV_v2.eps',format='eps',bbox_inches='tight')
plt.savefig('25meV_v2.pdf',format='pdf',bbox_inches='tight')
plt.show()
